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Transverse self-interactions within a "line" bunch of electrons moving in an arc of a circle have 
been studied extensively by us from an electrodynamical viewpoint. Here we treat an electron bunch 
with a given vertical (i.e. perpendicular to the orbital plane) size by a generalization of our previous 
work to the case of a test particle with vertical displacement interacting with a line bunch. In fact, 
since a bunch with vertical extent can always be thought of as a superposition of displaced charge 
lines, all the relevant physical aspects of the problem are included in the study of that simple model. 
Our generalization results in a physically meaningful and quantitative explanation of the simulations 
obtained with the code TraFiC 4 as well as in successful cross-checking of the code. 

PACS numbers: 29.27.Bd, 41.60.-m, 41.75.Ht 



I. INTRODUCTION 

Nowadays the production of high-brightness electron 
bunches constitutes one of the most challenging and in- 
teresting activities for particle accelerator physicists. 

The growing interest of the community in this kind of 
beams is justified by important applications such as self- 
amplified spontaneous emission (SASE)-free-electron- 
lascrs operating in the x-ray region (see, amongst oth- 
ers, 0)- Similar bunches are also being considered for 
production of femtosecond radiation pulses by simpler 
schemes based on Cherenkov and Transition Radiation 
□ 

One of the challenges faced by physicists involved in 
this area is constituted by the presence of self-field in- 
duced collective effects. These effects may spoil the 
brightness of the electron beam. 

Self-fields obey the usual Lienard-Wiechert expres- 
sions; the fields generated whenever an electron bunch 
undergoes a motion under the influence of external forces 
are different from the usual space-charge self-interaction, 
and are usually neglected in cases in which one does not 
deal with high-brightness beams. 

The equations for the field form, together with the dy- 
namical equations for the particle motion, a formidable 
self-consistent problem. Its complete solution, that is the 
particle evolution, is obtained only when one is able to 
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solve simultaneously the equations for the fields (elec- 
trodynamical problem) and the equations of motion (dy- 
namical problem): this task can be performed with the 
help of self-consistent computer simulations. An example 
of such codes is given (see by the program TraFiC 4 , 
which will be used throughout this paper. 

TraFiC 4 has been used to predict self-field related ef- 
fects in the bunch compressor chicane to be used before 
the main linac for the XFEL at DESY and it is, at the 
time being, the only fully developed code used for XFEL 
modelling both at DESY and SLAC. Results show that 
the projected transverse emittance of the bunch grows 
from 0.8 mm mrad after the injection to the significantly 
enhanced value of 2.6 mm mrad after the compressor 

Transverse dynamics is addressed by the code in two 
steps: first, the transverse electromagnetic forces, which 
are well defined and measurable physical quantities, are 
calculated separately and, second, the equation of motion 
is solved in a self-consistent way. 

From this viewpoint a thorough understanding of 
transverse electromagnetic forces by themselves is of 
paramount importance since it allows to gain confidence 
in full simulation results. The behavior of these forces 
can be obtained numerically by using the electromag- 
netic solver of the code for a given source distribution 
evolving rigidly along a certain trajectory (i.e. a zero en- 
ergy spread is assumed throughout the bunch evolution) , 
but numerical results alone are not sufficient in order to 
reach a full understanding of how electromagnetic forces 
act. Consider for example the plot in Fig. obtained 
by using the electromagnetic solver of TraFiC 4 . 

The figure shows the normalized radial force (i.e. the 
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FIG. 1: Radial force normalized to e 2 X/(4:neoR), F, felt by a 
particle with different vertical displacements h positioned at 
the center of a line bunch 200 /ira long as the bunch enters a 
circle of radius R = 1 m. Here 7 = 100. These simulations 
were obtained by means of the code TraFiC 4 . 



force normalized to e 2 \/ (A-keqR) in the direction orthog- 
onal to the test particle velocity and still lying on the 
bending plane) felt by a test particle. Here A is the bunch 
density while h is the vertical (orthogonal to the bending 
plane) displacement with respect to the center of a one- 
dimensional bunch with rectangular density distribution 
function. 

The plot shows interesting characteristics: a very sharp 
feature at the injection position, a transient and a steady 
state, which we can observe but for which there is no intu- 
itive physical explanation. The electromagnetic problem 
is, of course, only a part of the formidable evolution prob- 
lem but its understanding is perse an important physical 
issue. 

Analytical investigations of this simplified, although 
very important situation are therefore needed in order to 
see the meaning of the simulation results, besides provid- 
ing, from a more practical viewpoint, important cross- 
checks. We will show that through this kind of study 
we can understand the features of Fig. ^ from a physi- 
cal viewpoint, deal with characteristic times and lengths 
and, as a very practical result, validate the simulation 
results by producing analytical cross-checks. 

The self-interaction in the longitudinal direction (par- 
allel, at any time, to the velocity vector by definition) 
is responsible for the energy exchange between the sys- 
tem and the acceleration field and for all CSR (Coherent 
Synchrotron Radiation)-related phenomena, which have 
been studied extensively elsewhere (see, amongst others, 

!... El). 

Self-forces in the transverse direction were first ad- 
dressed, in the case of a circular motion, and from an elec- 
trodynamical viewpoint, in (l4|. Further analysis f[IT|. 

0) considers, again, the case of circular motion 
both from an electrodynamical and a dynamical view- 



point. We proposed a fully electrodynamical analysis of 
a bunch moving in an arc of a circle in 0, where, for 
the reasons explained before, we were not interested in 
the full evolution problem. Nevertheless, in that paper 
we addressed the issue of understanding the electromag- 
netic interactions only in the framework of a ID model. 
Before our study in [lj there was no explanation for the 
simulation results shown in Fig. ^ The work in 
provided a qualitative explanation of these results: such 
qualitative character of the explanation is due to the fact 
that the model in is one-dimensional, and does not al- 
low any quantitative investigation when (as in Fig. ^| a 
vertical displacement is introduced. Still it was possible 
to demonstrate that the sharp feature near the injection 
point is due to interactions of the test particle with par- 
ticles in front of it, while the transient part is due to the 
interaction with particles behind it. 

In this paper we aim at an extension of the model pro- 
posed in 1] which allows a fully electrodynamical analysis 
of a bunch endowed with a vertical dimension. Since a 
bunch with vertical extent can be always thought of as 
a superposition of displaced charge lines, all the relevant 
physical aspects of the problem are included in the study 
of a simpler model, constituted by a one-dimensional 
line bunch and a test particle with vertical displacement, 
which is the situation studied in Fig. ^ I n this article, 
our explanation of the features in Fig. ^wiU be eventu- 
ally refined to a quantitative level. 

The paper is organized as follows. We first treat, in 
Section [HJ the transverse interaction between two par- 
ticles moving in an arc of a circle, supposing that one 
of the two particles has a vertical displacement h with 
respect to the source. By integration of our results we 
consider, in Section IlIII a stepped-profile electron bunch 
interacting with a test particle with vertical displacement 
entering an arc of a circle and we discuss all the charac- 
teristic lengths involved. Finally, in Section ITV1 we come 
to a summary of the results obtained and to conclusions. 



II. TRANSVERSE INTERACTION BETWEEN 
TWO ELECTRONS 



We will first address the case of two electrons moving 
in an arc of a circle of radius R, in such a way that one of 
the two is displaced, with respect to the reference trajec- 
tory, by a quantity h in the vertical (perpendicular to the 
bending plane) direction. We will see that, as concerns 
the interaction in the radial direction (in the bending 
plane) it does not matter which particle is endowed with 
this displacement. 

The electro-magnetic force which one of the two par- 
ticles (designated with "T", i.e. the test particle) feels, 
due to the interaction with the other one (designated 
with "S", i.e. the source particle), is given by 



F(r Tl t) = eE(r T ,t) + ec(3 T x B(r T ,t), (1) 
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where rx is the position of the test particle, e is the 
electron charge with its own (negative) sign, (3t is the 
velocity of the test particle normalized to the speed of 
light, c, while E(rT,i) and B(rT,t) are, respectively, 
the electric and the magnetic field generated at a given 
time t by the source particle S, at the position of the test 
particle T, namely 



E(r T ,t) 



n-As 



47re [is R 2 ST (1 - n ■ f3 s y 



n x 



(fi - /3 S ) x (3 S 



R ST (l-n-(3 s y 



(2) 



two points and 7s is the usual Lorentz factor referred to 
the source particle at the retarded time t'. 

In contrast to the case studied in the transverse 
(by definition, orthogonal to /3t) projection of F(rT,i) 
has now components both in the bending plane and per- 
pendicular to it. We will introduce, therefore, two unit 
vectors §h, in the direction perpendicular to the bend- 
ing plane, and e p , in the radial direction, i.e. orthogonal 
to /3t and lying in the bending plane. Of course, the 
transverse component of F (rx, t) can still be written, fol- 
lowing [lj , as the sum of contributions from the velocity 
("C", Coulomb) and the acceleration ("R", Radiation) 
fields, namely 



and 



F±(r T ,t) = F ±c (r T ,t) + F ±R (r T ,t), 



(4) 



B(r T ,t) = -fix E{r T ,t) . 



(3) 



where 



Here (3s and (3s are, respectively, the dimensionless ve- 
locity and its time derivative at the retarded time t' , 
Rst is the distance between the retarded position of the 
source particle and the present position of the test elec- 
tron, fi is a unit vector along the line connecting those 



F± c (r T ,t) 



and 



n ± (l-/3 s -/3 T )-/3 iS (l-fi-/3 T ) 
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F± R (r T ,t) 
(3 ±T (1 - fi • (3 T ) + i 



4tt£ 
/St -Pa 
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(6) 



A. Tail-Head interaction: case of two particles in 
circular motion 

We will first consider the situation in which the test 
particle is in front of the source. In this case, one can 
refer to Fig. [21 for all the possible configurations of the 
present position of the test electron and the retarded po- 
sition of the source with respect to the arc. The vertical 
displacement h of one of the two particles is to be imag- 
ined in the direction perpendicular to the figure plane. 



Let us start with the steady state case in Fig. 03. We 
can define with As the curvilinear distance between the 
present position of the test and of the source particle; (f> 
will indicate the angular distance between the retarded 
position of the source and the present position of the test 
electron, and it will be designated as the retarded angle. 
Finally, h will be the vertical displacement. 

In the following we will assume (3s — Pt — (3- There- 
fore we can write Eq. (J5J and Eq. © in the following 
way: 



_LC 



eh [h(l — (3 2 cost 



[2i?sin 2 (</>/2)(l + (3 2 ) - /Jsin^/i 2 + AR 2 sin 2 (0/2)) 1 /2] 



47T£o7 



[{h 2 +AR 2 sm 2 ((f>/2)) 1 / 2 - 2/3i?sin(</>/2)cos(</>/2)]' 



(7) 
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FIG. 2: Relative configuration of the retarded source point 
S' and the test point T for a system of two electron passing 
a bending magnet. The test particle is understood to be dis- 
placed by a quantity h in the vertical direction. The figure 
shows only the projection of the motion on the bending plane. 



e 2 /3 2 c 

47T£ 

j a /o ^ e h [fr(l - P 2 cos0)] + 2e p sin(<ft/2) [flsin(0/2)(l + 1 ) - /3cos(^/2)(fe 2 + 4fl 2 sin 2 (</»/2)) 1 / 2 ] 
[ [(/i 2 + 4i? 2 sin 2 (<^/2)) 1 /2 - 2/3i?sin(^/2) cos(0/2)] 

e h [ft/3 sin 0] +e p sin - cos0(ft 2 + 4i? 2 sin 2 ^)) 1 / 2 ] 1 
R [(ft 2 + 4i? 2 sin 2 ( ( /)/2)) 1 /2 _ 2/3i?sin(0/2) cos(0/2)] 2 J ' 



(8) 



The retardation condition linking As, ft and <f> reads 



As = Rcp-/3 



h l + 4R 2 sin 2 



1/2 



(9) 



As one can readily see by inspecting Eq. @ , when we im- 
pose reasonable values for As <C R and ft <C R we obtain 
corresponding values of <fi <C 1 • We will therefore assume 
<C 1 throughout this paper, and verify a posteriori 
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the validity of this assumption when studying particular 
situations. Note that, by fixing < 1 we keep open 
the possibility of comparing (j> with the synchrotron ra- 
diation formation angle 1 /7 (note that a deflection angle 
smaller or larger than 1/7 is characteristic of the cases, 
respectively, of undulator or synchrotron radiation) . We 
can therefore expand Eq. Q and Eq. © to the second 
non-vanishing order in <f> thus obtaining 



4tte R 2 



(10) 



and 



e 2 7 3 
Ane R 2 

where we define 3>c and $r as 



(11) 



4ev, 



2 -4ft 2 /^ 



l + 2 /4 + /i 2 /0 2 



(12) 



and 



The latter can be written down in dimensionless form too 
as 



As = 

2 24 2 



(15) 



which explains the choice of the definition for h. 

The reader may easily check that Eq. lfT2jl. Eq. i(T3|) 
and Eq. lfT3|l (as well as Eq. J7J) and Eq. ©) are gener- 
alizations of the expressions given in [l| by taking their 
limit for h — ► 0. 

The following expression, which is valid for the total 
transverse force felt by the test particle can be then triv- 
ially derived 



e 2 7 3 
4vre i? 2 



(16) 



where $ is defined by 



*(0) = $r(0) + *c(</>), 



B. Tail-Head interaction: case (a) 



(17) 



*r(<£) 



1 - 4/i 2 /<^ 



2 - 4> 2 - 3h 2 + Ah 2 /(j> 2 



i + 4> 2 /4 + h 2 /<f> 2 



(13) 



Here and above 
ready treated in 



= ycj). This normalization choice, al- 
, is quite natural, I/7 being the syn- 
chrotron radiation formation-angle at the critical wave- 
length. In the derivation of Eq. I|1U|) and Eq. <|11[) (and in 
the following, too) we understood <f> ^ I/7, which is jus- 
tified by the ultrarelativistic approximation. Moreover 
we defined h = h^ 2 /R; in order to understand this defi- 
nition we first write down the retardation condition Eq. 

in approximate form. Since we are already working 
in the limit As/R <C 1 and <§; 1, then in order to have 
As > 0, h < R(j) must be satisfied. We will automatically 
recover from our results that the assumption h <C Rip is 
sufficiently good for our purposes. Therefore Eq. @ can 
be approximated with: 



As = (1 - /3)i?</> 



R<j) 3 f3h 2 
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2R<j) 



(14) 



Let us now consider the other cases depicted in Fig. 
J2Kj c and d. While the case in Fig. 03 deals with the 
steady state situation in which the present position of the 
test and the retarded position of the source are both in 
the bend, Fig. c and d deal with transient situations 
in which we can find the retarded source and the present 
test in the straight line before and after the bend too. 

Consider the situation in Fig. In this case, under 
the assumption h <C (y + R<fi) > the retardation condition 
reads 



2 



^/3(y + ^/A) 
y + cj> 



(18) 



where we introduced the normalized quantity y = yj/R, 
which is just y/R normalized to the synchrotron radiation 
formation angle, I/7. 

In this case the source particle is only responsible for 
a velocity field contribution, therefore F± = F±c- By 
direct use of Eq. (JSJ, one can find the exact expression 
for F ± 



e 2 e h {h[l- P 2 cos (/)}} +e p {-R ST P sin <j)+[R(l + (3 2 ){1- cost/)) + y sin0] } 
± ~ 4tt£o7 2 {Rst - (3{y + i?sin0)} 3 

I 

Rst being given by Expanding the trigonometric functions in Eq. (|19|) and 



Rst = [{y + Rsmc/)) 2 + R 2 (l -coscp) 2 + h 2 ] 1/2 . (20) 
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using normalized quantities one finds: 



1.2- 




0.00 0.02 0.04 0.06 0.08 0.10 



Position after injection (m) 



FIG. 3: Normalized radial force F for a two-particle system 
entering a hard-edge bending magnet as a function of the 
position after injection. Results for different values of h are 
shown. The solid lines show analytical results; the circles 
describe the outcome from TRAFIC 4 . Here As = 5.0. 



e 2 4-7 3 

47T£ R 2 



(y + 4>y 



|e h 




+ e p 


y 2 <p + y{cf> 2 + <^ 4 /2) + 5 /4 - h 2 <f) 


} 




(y + 0) 2 + 4 /4 + h 2 
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(21) 



It can be easily verified that, as it must be, Eq. H21|) 
reduces to Eq. I jlOjl in the limit y — ► 0. Also, the reader 
may check that in the limit h — > 0, Eq. Ij21|l reduces to 
already derived expressions in 0. 

Let us now define the normalized radial force F = 
F p /[e 2 /(ATTe RAs)). It is possible, by means of Eq. (|2~TJl . 
to plot F as a function of the position after the injection 
(defined by the entrance of the test particle in the hard- 
edge magnet) for a fixed value of As = Asj 3 /R = 5.0 
and different values of h. In Fig. |3we compare such a 
lot with numerical results from the code TRAFIC 4 (see 



Note that, as already pointed out in [jj], at the posi- 
tion which corresponds to the entrance of the retarded 
source in the magnet there is a discontinuity in the plots. 
This is linked to our model choice, and it is due to the 
abrupt (hard edge magnet) switching on of the accelera- 
tion fields. 

As general remark to Fig. [3| and, in fact to Fig. 0] 
Eland Fig^l the perfect agreement between our calcu- 
lations and numerical results by TRAFIC 4 provides, per 
se, an excellent cross-check between our analytical re- 
sults (with their assumptions and applicability region) 
and simulations. 



C. Tail-Head interaction: case (c) 

We will now move to the case depicted in Fig. [2J;, in 
which the source particle has its retarded position inside 
the bend and the test particle has its present position 
in the straight line following the magnet. We will de- 
fine with x the distance, along the straight line after the 
magnet, between the end of the bend and the present po- 
sition of the test particle. Here <fi m = 7</> m , 4>m being the 
angular extension of magnet, and x = 'yx/R, the reason 
for this normalization choice for x being identical to that 
for y. 

The retardation condition reads 



As 



if) + x ip 3 if; + Ax 
2 + 24 ^ + £ 



2(if) 



(22) 



In contrast with the case of Fig. here we have contri- 
butions from both velocity and acceleration field. Again, 
by direct use of Eq. JSJ and Eq. JBJl one can find the 
exact expression for the transverse electromagnetic force 
exerted by the source particle on the test particle 



F± = F±c + -Fj_r , 



(23) 
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where 



e 2 e h [(1 - /3 2 cosip )h] + e p [i?(l -cosip)(l - (3 2 cosip) - (3smip(R S T ~ Px - (3Rsm ip)] 
47T£o7 2 [Rst - /3 x cos ip — (3R sin ip] 



(24) 



and 



^ §h [/5 2 /i(l — /3 2 cos ^>)(i2 + xsin ^ — i?cosi/')] 

±R AtteqR I [i? ST - fa cos V - /3-R sin ijj} 3 

e p f3 2 [(R + xsimp - Rcosip)(R(l - cosi/>)(l - /3 2 cos?/>) - /3 sin ip{Rsr ~ fix + /3i? sin ?/>))] 



+ 



[Rst — fix cos ip — (3R sin ip] 3 
§h/3 2 [/isin?/)] + e p [R(l — cosip)f3 2 s'mip — (3 cosip(RsT — fix — ptRsmtp) 



[Rst — fix cos ip — (3R sin ip] 
I 



(25) 
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Rst being now 



R ST = ((x + RsiniP) 2 + R 2 {1 ~ coaiP) 2 + h 2 ) 1/2 (26) 



FIG. 4: Normalized radial force F for a two-particle sys- 
tem leaving a hard-edge bending magnet as a function of the 
position after the ejection. Results for different values of h 
are shown. The solid lines show analytical results; the circles 
describe the outcome from TRAFIC 4 . Here As = 5.0. 



Expanding the trigonometric functions in Eq. (|24|) and 
Eq. (|25|l . and using normalized quantities one finds: 



2e 2 7 3 2e h h{2 + tp 2 )(x + </>) 3 + e p V3(i + ip) 2 -2x 2 + x (tP 3 - 2tpj + ip 4 /2 - 2fi 



4ire Q R 2 



(x + iP) 2 + {iP 2 /4)(2x + iP) 2 + h 2 



(27) 



and 



2e 2 -v 3 - f- 2g h 

f ^=^r^ + M — 



hip(x + ip) + e p x 2 + xip(2 - ip 2 ) + ip 2 - (3/4)V> + h 2 



{x + ip) 2 + (^ 2 /4)(2i + ^) 2 + h 2 
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2e h (i + ?A) 2 (2 + ^ 2 ){x + 4>/2)4>h + e p (x + tp)(x + ?/>/2)^ 2 [ - 2h 2 - 2x 2 + xtp{-2 + ip 2 ) + tp 4 /2] 



(x + </>) 2 + (^ 2 /4)(2i + vo 2 + h 2 



(28) 
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FIG. 5: Time variation of a synchrotron radiation pulse gen- 
erated by a highly relativistic electron moving in a circle as 
seen by an observer in the orbital plane 



Similarly to the latter case, it can be easily verified that 
Eq. J23 and Eq. ggj reduce to Eq. (JTUJ and Eq. (jTT|) . 

respectively, in the limit x — > 0. Moreover, in the limit 
/i — * 0, they reduce to already known expressions given in 
. Again, one can plot the normalized radial force F as a 
function of the position after the ejection, defined by the 
exit of the test particle from the hard-edge magnet, for 
different values of h and a fixed value of As = Asj 3 / R = 
5.0. In Fig. 01 we compare such a plot with numerical 
results from TRAFIC 4 . 

At the position which corresponds to the exit of the 
retarded source from the magnet there is a discontinuity 



in the plots. This, again, is linked to our model choice, 
and it is due to the fact that, for particles on axis, af- 
ter the retarded source has left the magnet there is only 
Coulomb repulsion along the longitudinal direction. 

It is suggestive to notice the resemblance of the peaks 
shown in Fig. 21 with half of the time pulse of the ra- 
dial electric field from usual synchrotron radiation pro- 
cess (see Fig. and 0). This is not a coincidence. The 
test particle is, indeed, far away from the source with re- 
spect to the formation length -R/7 3 and the magnetic field 
contribution to the Lorentz force can be expected to have 
the same behavior of the electric field contribution, since 
B = (n x E)/c. The only difference is that the observer 
is now "running away" from the electromagnetic signal 
which will result in a spreading of the pulse of about a 
factor (1 — P)^ 1 . Since R — 1 to and 7 = 100 we expect 
the pulse to be long about (i?/7 3 )(l - ~ 10~ 2 to. 
D. Tail-Head interaction: case (d) 

The last Tail-Head case left to discuss is depicted in 
Fig. the source particle has its retarded position in 
the straight line before the bend, and the test particle 
has its present position in the straight line following the 
magnet. The retardation condition reads 



+ x + y 



- 4(.t + y) + 12xy/<p m - h 2 /2 



JJ 



.(29) 



In this case we have only velocity contributions. The 
exact expression for the electromagnetic transverse force 
on the test particle is 



e 2 £?h [sin 0(1 — (3 2 cos</> m )] + e p [sin 5 cos#(l — f} 2 coS(f) m ) — /3sin0 m (l — cosfScos^)] 
F± — — ; — (30) 



47T£o7 



R ST [1 — [3 sin 4> m sin S cos 9 — [3 cos <j> m cos S cos 9]' 



where 



smf 



Rst 



(31) 



sin S = 



R + y sin (f> m — R cos < 
Rst cos 9 



(33) 



cos 9 



{R 2 ST -h 2 f' 2 
Rst 



and 



(32) 



cos S 



x + y cos (j) m + R sin < 
Rst cos 9 



(34) 
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where Rst can be retrieved by the latter two equations finds the following approximated expression for F±: 
and some trivial trigonometry. 

Once again, expanding the trigonometric functions in 
Eqs. I|o0[l ... (|33[) and using normalized quantities one 



F± ~ - 6 D2 47 3 (f + y + 4> m f$ m x 



eh 


(2/<f> m + 4> m ){x + y + 4> m )h 


+ e p 


-x 2 + y 2 + {4>l)xy + 4(&/2 - m ) + y($j2 + r> m ) + &/4 ~ h 2 




£(l + &)+0 + &n + &/3 


- (^712) [\2xy + 4(£ + y)0 m + &] + /j 2 }' 



• (35) 



It is easy to verify that Eq. i|o5|l reduces, respectively, 
to the steady state (Eq. ljl"fj|) 1 when a; = and y = 0, 
to the transient case in Fig. |2K when a; = (Eq. I|21|)) 
and to the transient case in Fig. when y = (Eq. 
iJUJl). Moreover the reader may verify that, in the limit 
h — > 0, Eq. (|35fl reduces to already known results in 
1] . Following the treatment of the transient situations 
in Fig. and in Fig. [2t we plot, for this case too, 
a normalized expression for the transient force, i.e. the 
usual F, as a function of the curvilinear position of the 
test particle (s = indicates the entrance of the magnet) 
for different values of h and a fixed value of As = Asj 3 / R 
and magnet length. In Fig. we compare our analytical 
results with numerical results from TRAFIC 4 , for a fixed 
value of As = 5.0 and <j) m = 1 and for different values of 
the vertical displacement h. 
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E. Head-Tail interaction 

Finally, we can deal with the situation in which the 
source particle is ahead of the test electron, i.e. As < 0; 
we will talk, in this case, of head-tail interaction. On 
the one hand it is evident that, when As < 0, the source 
particle is ahead of the test electron at any time; on the 
other hand it is not true that the present position of the 
test particle is, in general, ahead of the retarded position 
of the test particle. As already suggested in 0, if As < 
and, approximatively, |As| < h the test particle overtakes 
the retarded position of the source before the electromag- 
netic signal reaches it. In this case, although we may still 
talk about head-tail interaction, since As < 0, its real 
character is very much similar to the case As > 0, just 
treated in the previous subsections, in which the electro- 
magnetic signal has to catch up with the test particle. In 
order to understand the physics involved in this situation 
it is sufficient to study the cases As < and |As| > h. 
The case As < with |As| < h will be treated from 
a qualitative viewpoint only: its quantitative analysis, 
which may be interesting to perform for the sake of com- 
pleteness, presents stronger mathematical difficulties and 



FIG. 6: Normalized radial force F for a two-particle system 
crossing a hard-edge bending magnet as a function of the 
position of the test particle inside the magnet in the case 
of a short magnet r/> m <C 1. Results are shown for different 
values of h. The solid lines show analytical results; the circles 
describe the outcome from TRAFIC 4 . Here the normalized 
distance between the two particles is As = 5, while 4>m = 1. 



is left for future study. Anyway a qualitative treatment 
of this situation is enough to reach a full understanding 
of the interaction physics, which is our goal here. 

Consider therefore the case As < with | As| > h. It is 
important to note that, once the steady state case is stud- 
ied, the situation in which the source particle is ahead of 
the test electron can be treated immediately for all three 
(a, c and d) transient cases in Fig. [5] (of course, with 
respect to the figure, test and source particle exchange 
roles) on the basis of the steady state case alone. In fact 
in that situation, as it has already been said in we can 
assume the retarded angle </> small enough (the test par- 
ticle "runs against" the electromagnetic signal) that the 
actual trajectory followed by the particles is not essential 
and one can use the steady state expression to describe 
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also the transient cases. It can be shown that the only 
important contribution from the source particle comes 
from the acceleration part of the Lienard-Wiechert fields. 
Within our approximations, the only non-negligible con- 
tribution is present in the situation (again, with the roles 
of test and source particle inverted) depicted in Fig. [2t 
and Fig. the latter just being the steady state case. 

Let us deal with the steady state case of head-tail in- 
teraction in the case \As\ > h. As already discussed in 
0, the difference with respect to the situation in which 
the test electron is in front of the source is that the test 
electron "runs against" the electromagnetic signal emit- 
ted by the source, while in the other case it just "runs 
away" from it. Therefore the relative velocity between 
the signal and the test electron is equal to (1 + (3)c, in- 
stead of (1 — j3)c as in the other situation. Hence the 
retardation condition reads 



As~R<j) + p [R 2 



I 1/2 



(36) 



or, solved for 



As 



As 2 I 1 



r 



1/2 



(37) 



Note that, in the asymptotic for As ^ 7/1 we recover the 
result in [l| for the one-dimensional case. 

In the general case, /3s is almost parallel to (and equal 
to) (3t and antiparallel to the projection of n on the 
bending plane: it turns out that the only important con- 
tribution to the radial force on the orbital plane is given 
by the second term on the right side of Eq. © , and it is 
easy to check that 



Fi 



(38) 



AneoRAs 

It may be worthwhile to underline that the force in Eq. 



(|38|l is, structurally, identical to the result in pj although, 
of course, As is now a function of h . 



III. TRANSVERSE INTERACTION BETWEEN 
AN ELECTRON AND A BUNCH ENTERING A 
BEND FROM A STRAIGHT PATH 

In the previous Section we dealt with all the possible 
configurations for a two-particle system moving in an arc 
of a circle. Now we are ready to provide a quantitative 
explanation of Fig. 

After the discussion, in Section [Hj about head-tail in- 
teractions within a system with the test particle behind 
the source electron, one is led to the qualitative conclu- 
sion that, within an electron bunch, interactions between 
sources in front of the test particle and the test particle 
itself are important and, in general, they must be re- 
sponsible, at the entrance and at the exit of the bending 
magnet, for sharp changes in the transverse forces acting 
on the test electron. The quantitative change depends, 
of course, on the position of the test particle inside the 
bunch: the extreme cases are for the test particle at the 
head of the bunch, where there are just interactions with 
electrons behind the test particle (no head-tail interac- 
tions), and for the test particle before the tail of the 
bunch, at a distance |As| > h, where all the sources are 
in front of it (only head-tail interactions in the regime 
|As| > h). 



A. Head-Tail interaction 

Consider now the situation in which a one-dimensional 
line bunch interacts with a particle positioned at the cen- 
ter of the line but vertically displaced by a quantity h. 
As already said we will discuss, analytically, the head- 
tail interaction for |As| > h only. Since the trajectory 
followed by the bunch is not important we can simply 
integrate Eq. I|38ll and find the contribution 



J 



p HT 



(0) ~ <^ e 2 Ao/(47r£ 0j R) ln[As max /(As 
I e 2 A /(47T£ 0j R) ln[As max //i] 



i?0(l + /?))] 



cf>< 
< <j> < $(/i, As max ) 



(39) 



r 



where "HT" stands for "head-tail" and As) is the 
solution, in 0, of Eq. (|36() at vertical displacement h and 
longitudinal distance As. 

When a bunch longer than the vertical displacement 
h enters the magnet, the particles in front of the bunch 
will interact with the test electron following Eq. (1391) . 
which models the radial interaction on the basis of Eq. 
(|38p. Nevertheless, Eq. I|38|l cannot describe the situ- 



ation when the sources at a distance shorter than the 
vertical displacement h begin to interact with the test 
electron. As a result the head-tail interaction has two 
characteristic formation lengths. The first one indicates 
the distance that the test electron travels from the mo- 
ment it is reached by the electromagnetic signal from 
the first particle entering the bend, till the moment it 
is reached by the electromagnetic signal emitted as the 
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FIG. 7: Normalized radial force F from the head of the bunch 
(As < 0) to the tail as the bunch progresses inside the bend 
(z = corresponds to the injection of the test particle in the 
magnet). Here R = lm, 7 = 100, h = 10/im and the bunch 
length is 200/im.The test particle is located at As = 0. 

particle at As = —h enters the bend. This is given by 
L-y = /3cR((f>2 — 4>\) — (3cR(j)2, 4>2 an d <j>i being the so- 
lution of Eq. (|37l) when As = lb/2 and As = /3h, re- 
spectively (the reader will recognize that <j>\ = 0). In 
the limit for h « lb/ 2 Eq. J37J), substituted in the ex- 
pression for L±, gives the rule of thumb L\ ~ Z&/4. The 
second characteristic length is given by the distance that 
the test electron travels from the moment it is reached by 
the electromagnetic signal from the particle at As = —h 
till the moment it is reached by the electromagnetic sig- 
nal emitted by the particle with As = as it enters the 
bend. This can be estimated to be, roughly, equal to 
7/1, at least when h is not too large. In fact, in order 
to know the present angular position of the test particle 
when As = one should solve the retardation condition 
in (f> 

R(/> = /3[h 2 + sin 2 (0/2)] 1/2 (40) 

In the limit in which sin(</>/2) ~ <j>/2 we have 4> ~ jh/R 
and L2 — 7/1. 

This reasoning can explain pretty well the features in 
Fig. [7| The radial force grows exactly as described in Eq. 
(|39H starting from about z — —lb/4: — 50 /im until z — 
(which is roughly the first formation length). Note that 
with the source located at As = h corresponds a retarded 
angle = 0, which means that the electromagnetic signal 
from this source, which is emitted at z — 0, will reach the 
test particle when this has also the position z = 0. After 
z = 0, the test particle begins to interact with the sources 
located at As > — h while the ones at As < —h reach a 
steady state, in the sense that they keep on interacting 
in the same way with the test electron. 



FIG. 8: The same as Fig. |7|for a different range. Normalized 
radial force from the head of the bunch (As < 0) to the tail 
as the bunch progresses inside the bend (z — corresponds 
to the injection of the test particle in the magnet). Here 
R — lm, 7 = 100, h — Wfim and the bunch length is 200/im. 
The test particle is located at As = 0. Curve A corresponds to 
the simulation by TraFiC 4 , curve B is the contribution, given 
by Eq. 1391 of the particles with As < —h; A-B is the difference 
between the two curves, ascribed to the contribution from the 
particles with As > —h. 



This is illustrated by Fig. |H1 where the normalized 
radial force F is plotted. 

The difference between the simulation results and the 
analytical estimation in Eq. Ij39(l gives, quantitatively, 
the radial interaction due to the sources at As > —h. 
For h — 10 fim and 7 = 100, as in our case, one expects 
a second formation length equal to 10 -3 m which is ex- 
actly what one gets: the actual data show, in fact, that 
a maximum is reached when z = 0.001 m at F ~ 7.098. 

Fig. [S]also shows the steady state regime, when all the 
particles in front of the test particle interact with the lat- 
ter from inside the bend. The solid line in Fig. |S1 shows, 
once again, the interaction with the sources located at 
As < —h. Finally, in Fig. Elwe present the results from 
TraFiC 4 for several values of h. Of course, when h> h, 
i.e. h > 100 fim there are no sources characterized by 
As < —h at all. 

One may expect that the system enters the steady state 
at z ~ jh, which is correct for h = 1 jim, h = 10 /im 
and, by figure inspection, for h = 100 fim. Nevertheless, 
in the case h = 1 mm and h — 10 mm the system en- 
ters the steady state at, respectively, z ~ 0.054 m and 
z ~ 0.186 rn, which are clearly smaller than 7/1. The rea- 
son for this apparent discrepancy is due to the approxi- 
mation sin(0/2) ~ <j)/2 which has been used to derive the 
rule of thumb L 2 — "fh starting from the retardation con- 
dition Eq. lEOJ. A comparison between the rule of thumb 
proposed before (dashed line) and the real solution of the 
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FIG. 9: Normalized radial force from the head of the bunch 
(As < 0) to the tail as the bunch progresses inside the bend 
(z = corresponds to the injection of the test particle in the 
magnet). Here R — lm and 7 = 100. Here we plot the results 
from TraFiC 4 for several values of h. The bunch length is 
200/im. The test particle is located at As = 0. 

retardation condition (solid line) is given in Fig. 1101 one 
can easily see that, when h = 1 mm, (j> ~ 0.054. In the 
same way, at h = 10 mm, <f> ~ 0.186 (remember that 
R = lm). 

By comparison, in Fig. [7| between the contributions 
to the head-tail interaction in the cases As < —h and 
—h< As < the reader can see that these are compa- 
rable in magnitude, although h is much shorter than the 
bunch length. A qualitative explanation of this fact can 
be given. When As <C <f>R, so that it can in fact be ne- 
glected in the retardation condition, one gets <fi ~ 7/1/i?. 
Now, the force felt by the test particle in the zeroth order 
in 4> is proportional to 7 3 /(i? 2 r/<) ~ j/(Rh): simple con- 



siderations suggest that this order of magnitude will not 
change appreciably up to As of order — /1/7. As a result, 
the region (0,-/1/7) gives a contribution, after integra- 
tion, proportional, approximately, to X/R. On the other 
hand, the maximum interaction in the region As < —h 
is of the order of e/(Rh), which also corresponds, after 
integration, to a magnitude proportional to X/R. This 
explains the reason why a small region —h< As < can 
give contributions of the same order of magnitude as the 
much longer interval in which As < —h. 



B. Tail-Head interaction 

We will now analyze the tail-head part of the interac- 
tion when As > 0. In the case where the bunch enters the 
bend we have contributions from retarded sources both 
in the bend and in the straight line before the bend. The 
contribution from the retarded sources in the magnet can 
4> 
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FIG. 10: Solution of the retardation condition (solid line) 
at As = as a function of h and comparison with <f> = 7/1 
(dashed line) 



be obtained by simple integration of Eq. l|lo|. an d reads 



lm 



e 2 X Q 
iireoR 



In 



A(2h 2 + t&J) 



A{2h 2 + 4> 2 min ) 



Ah 2 



b 2 

max 



b 4 

max 



4^ 2 + 4^2. + 



(41) 



where "m" is a reminder that the contributions treated 
by Eq. (|4*T|l are all from the " magnet" . All that is left to 
do now, is to investigate the values which </f> m i n and </> max 
assume. 

Let us first define with r) the normalized angular posi- 
tion of the test particle inside the bending magnet. Now 
define with (f>* the solution of the retardation equation 
As min = $*/2 + 0*724 - h 2 /{2$*). If (j>* < fj, the re- 
tarded position of the first source particle is in the bend- 
ing magnet, and m i n = 4>* ■ On the other hand, when 



(jf > 4> there are no contributions to the radial force from 
the bend. 

Next, we define with </>** the solution of As max = 
0**/2 + 0**724 - h 2 /{2cj)**) (in our case As max will be 
equal to one half of the normalized bunch length). Sup- 
posing cf>* < ij, if (/)** < r} too, then all the particles 
contribute from the bend, and max = <j)** . On the other 
hand, when </>** > f], we have a mixed situation, in which 
part of the particles contribute from the bend and others 
from the straight line before the magnet. In this case 
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FIG. 11: Normalized radial force acting on a test particle 
from a bunch with rectangular density distribution entering a 
hard-edge bending magnet as a function of the position of the 
test particle inside the magnet. The solid lines show analytical 
results; the circles describe the outcome from TRAFIC 4 . We 
chose As max = 100/im, 7 = 100, R = 1 m; graphs are plotted 
for several values of the parameter h. 

The contribution from the retarded sources in the 
straight path before the bend is given by 



F 



As 



R 

/y3 



F ± (y(As,h),<f>)dAs , (42) 



where "s" stands for "straight path", and where the ex- 
pression for F± in the integrand is given by Eq. (|21[l . It 
is convenient, as done before, to switch the integration 
variable from As to y. The Jacobian of the transforma- 
tion is given by (see 0) 



dAs _ (0 + j/) 2 + 4 /4 + fr 2 
dy ~ 2(4> + y) 2 



(43) 



After substitution of Eq. gSJ and Eq. in Eq. l|32)l. 

one can easily carry out the integration, thus getting 



F 



4> Uvmin + 20 + : 



Uvmax. + 20 + 3 

(4y2 nax + W) + 8y max + 402 + 4 



(4.4) 



As done before for 
the values of y mni and y 



m i n and max , we can now investigate 



Let us start with y min . First, we define with y* the so- 
lution of the retardation condition As m i n = (<fi + y*)/2 + 

(0 3 /24)(4y* + 4>)/(y* + 4>) - h 2 /(2y* + 20). If y* > 0, 
the retarded position of the first source particle is in the 
straight line before the bending magnet, and y mU i — y* ■ 
On the other hand, when y* < 0, the retarded position 
of the first source particle is in the bend, and y mm = 0. 

Next, we define with y** the solution of As max — (0 + 
y**)/2 + (0 3 /24)(4y**; + + 4 _ h 2 /{2y** + 20)) 

(again, in our case As max is just half the bunch length 
normalized to i?/7 3 ). Consider the case y** < 0: all the 
particles contribute from the bend, that is we entered the 
steady-state situation. In this case j/ max = ymm = 0. On 
the other hand, when y** > 0, we have again a mixed 
situation, in which part of the particles contribute from 
the bend and others from the straight line before the 
magnet. In this case y max = y** ■ 

The following step is to actually plot the radial force 
exerted on an electron by a bunch with rectangular dis- 
tribution entering a long bend. Our results, compared, 
once again, with simulations by TRAFIC 4 , are shown in 
Fig. ^2f° r a bunch length of 100 /xm, 7 = 100, R = 1 m 
and for different values of As. A perfect agreement is ob- 
tained with the results by TraFiC 4 . A first characteristic 
length is obviously given by the solution of the retarda- 
tion condition with As = h/2, and gives the position at 
which the test particle begins to feel a steady interaction 
from the tail sources. There is a second characteristic 
length more difficult to sec, though: a careful inspection 
of Fig. 1111 shows, in fact, a small irregularity (actually 
a discontinuity in the first derivative) in the curve for 
h = 1 at the position z ~ 1 cm. Such irregularities are 
present in all the curves in Fig. ^2 although one has to 
look carefully for them by magnifying parts of the plots, 
as shown in Figs. I12IA-D. It is this fact which actually 
suggests the presence of a second formation length. 

When the value of h is small the discontinuity is lo- 
cated, approximately, at z ~ 7/1: for example when 
h = 1 //TO, in Fig. 1121 ^. we have a discontinuity in 
the first derivative of the curve at z ~ 0.1 mm. Nev- 
ertheless, this value changes as we increase h. In fact, 
when h = 1 mm, in Fig. I12C . the discontinuity is at 
z ~ 5.4 cm, while when h = 1 cm, in Fig. I12L ). we find 
a value of z ~ 18.6 cm: these are the same numerical 
values found when discussing the entrance in the steady 
state of the Head- Tail interaction in the previous subsec- 
tion. The reader will remember that these are, in fact, 
the solutions of Eq. I|40|) . From a physical viewpoint, 
the solution of Eq. H40fl is the position at which the test 
particle begins to feel the electromagnetic signal from the 
source at As = 0. 

Before that point, the test particle feels interaction 
from particles behind it but only due to velocity fields, 
since the retarded positions of all the electrons behind 
the test one are not yet in the bend. After that particu- 
lar point, the force on our test electrons has a component 
due to the acceleration field too. This suggests that the 
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FIG. 12: Normalized radial force (analytical result) acting on a test particle from a bunch with rectangular density distribution 
entering a hard-edge bending magnet as a function of the position of the test particle inside the magnet. We chose As ma x = 
100/im, 7 = 100, R = 1 m; several curves for different values of h are shown. A clear discontinuity in the first derivative is 
visible. 



physical meaning of the presence of this second formation 
length is simply the switching on of the contribution of 
the acceleration field. As a last remark it is interesting 
to note that Tail-Head and Head- Tail interactions have 
a characteristic length in common, but for completely 
different reasons. 



IV. SUMMARY AND CONCLUSIONS 

In this paper wc presented a fully electrodynamical 
study of transverse self-forces within an electron bunch 
moving in an arc of a circle in the case the test particle is 
endowed with a vertical displacement h. We strived for 
a generalization of the results obtained in in order to 
obtain a better qualitative and quantitative explanation 
of the physics involved in the problem and, in particular, 
to explain the behavior of the self-interaction depicted in 



Fig.IH 

First we generalized the results in [lj in the case of 
a two-particle system. Then, by integration of these re- 
sults, the case of a line bunch and a test particle with a 
vertical displacement was studied. This case includes all 
the relevant physics present in the situation of a bunch 
with vertical size. 

Besides allowing one to generalize results obtained in 
, our study aimed at a physical understanding of the 
results by TraFiC 4 : we found that the bunch can be di- 
vided into four separate regions (over which one can inte- 
grate the two-particle interaction) corresponding to four 
different types of source-test interaction, namely head- 
tail with |As| < h, head-tail with |As| > h, tail-head 
with contributions from velocity fields alone, and tail- 
head with contributions from both acceleration and ve- 
locity field. We found that these regions correspond to 
four characteristic formation lengths, which can be deter- 
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mined quantitatively by simple analytical calculations. 

For the first, the third, and the fourth region, we could 
use relatively simple analytical results in order to de- 
scribe the situation and easily perform crosschecking with 
TraFiC 4 . The perfect agreement we found gives us much 
information: this constitutes, in fact, a reliable cross- 
check which provides, at the same time, very strong indi- 
cation that the code calculates the self-interaction in the 
proper way and that the approximations which we made 
in our analytical theory are correct. Moreover, we have 
physical explanations of the self-interactions in terms of 
formation-lengths and type of source-test interaction. 

Because of mathematical difficulties linked with the 
structure of the retardation condition we left the quan- 



titative discussion of the second region (— h < As < 0) 
for future work. Nevertheless we were able to understand 
the physical meaning of this region and to determine its 
formation length. 
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